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Abstract 

We derive Hamiltonian flow equations giving the evolution of the Lipkin 
Hamiltonian to a diagonal form using continuous unitary transformations. 
To close the system of flow equations, we present two different schemes. First 
we linearize an operator with three pairs of creation and destruction operators 
by reducing it to the z component of the quasi spin. We obtain the well known 
RPA-result in the limit of large particle number. In the second scheme we 
introduce a new operator which improves the resulting spectrum considerably 
especially for few particles. 
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Recently, a new method Q has been developed in condensed matter physics 
to find approximate solutions in many-body theory. This method is based on flow 
equations for the evolution of generalized Hamiltonian couplings with a fictitious 
parameter i running from zero to infinity. The evolution parameter I characterizes 
the interval of energies \Ei — Ej\ < 1 / y£, for which the Hamiltonian has sizeable 
non-diagonal matrix elements. In the limit i — > oo the Hamiltonian is diagonalized. 
The Hamiltonian flow equations have been successfully applied to the problems of 
electron-phonon coupling M , to BCS theory |3j and dissipative systems. It is natural 
also to try out this method on the Lipkin model which is a standard test case for 
approximation methods in nuclear physics. The Lipkin model is exactly solvable, 
since its Hamiltonian is composed of bilinear products of creation and destruction 
operators which have a quasi-spin representation. Recently, new interest in the 
Lipkin model has developed in the context of finite temperatures |j] and as a test 
of self-consistent RPA-type approximations ||. The Hamiltonian flow equations 
for the Lipkin model can be closed if one linearizes higher products of quasi-spin 
operators into a product of expectation values and a quasi-spin operator. When 
the number of particles N is very large we naturally reproduce the RPA results 
for N — > oo. Extra care is necessary to achieve satisfactory results for mesoscopic 
problems with iV = (8, 10.. 20) particles. We will give an example how to improve 
the simple linear scheme. 

In the Lipkin model N particles can distribute themselves on two levels, which 
are both iV-fold degenerate. The splitting of the levels is £o and the interaction Vo 
mixes states where two particles simultaneously move from the lower level to the 
higher level or vice versa: 

U = q£o Yl aa t a P° + 9^0 Y a tv a t'« a P'-° a P-° W 

Z o;P Z pp',cr 

In the Hartree-Fock ground state all particles are in the lower level, where they have 
the quasi-spin a = — 1. Therefore the total quasi-spin J z , 

J z = lY aa t a p^ ( 2 ) 

p,a 

has the value (J z ) = — |iV. 

The interaction term can also be expressed in terms of quasi-spin operators. The 
Hamiltonian 7i then has the form 

n = ^j z + v (j 2 + + j 2 _), (3) 

where ^ ^ 

J+ = —?k Y a p,+i a p -i> J- = ~7k Y a p~i a p+ti ( 4 ) 
yip \l i v 

and J z form an angular momentum algebra 

[Jz, J±] — ±^±; 

[J+,J_] = J z . (5) 
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Since [J 2 ,7i] = 0, the interaction mixes only states in the quasi-spin multiplet 
J = |iV. In the following we will consider a limited range of couplings such that the 
Hartree-Fock ground state remains stable, i.e. e = < 1, otherwise one would 
have to transform to a deformed basis || . 

The concepts of unitary transformations to diagonalize the Hamiltonian is well 
established. Negative energy states can be eliminated from the Dirac equation by 
unitary transformations |7], ||. The e s method || with a suitable ansatz for the 
operator S has been extremely successful to derive the hierarchy of two-, three- and 
more-body correlations in many-particle systems, especially nuclei. The new idea in 
the Hamiltonian flow equations is to make the unitary transformation infinitesimal, 
so that it can adjust itself best to diagonalize the Hamiltonian at each stage of the 
flow. 

In the field theory of critical phenomena it is advantageous to gradually inte- 
grate out momentum shells instead of eliminating the whole high momentum region 
in one step. In Hamiltonian theory the energy representation is the most useful rep- 
resentation for unitary transformations to eliminate non-diagonal matrix elements 
of shorter and shorter range. The resulting couplings in TC and therefore the Hamil- 
tonian itself will then carry an index £, labeling the stage of diagonalization. With 
f\ anti-hermitian the unitary transformed 7i has the following form 

U' = eme'V. (6) 

In the case of an infinitesimal fj = i]A£ we approximate the change 

Uil + M) - n(l) = [r)(t), H(£)]M (7) 

and obtain the Hamiltonian flow equation: 

dU{l) 



df 



W),W)\- 



In general one chooses rj proportional to the commutator of the diagonal part of 
the Hamiltonian 7Yd and n in order to decrease the magnitude of the off diagonal 
matrix elements during evolution: 

v (i) = [n D (£),n(£)}. (9) 

It is easy to follow the non-diagonal matrix elements of n in the course of evolution 
with L 

—Tf = ^(VikHkj ~ n ik r]kj) (10) 

and 

rjik = (Ei - E k )n ik (11) 

gives 

^ = -(E t - E 3 fU %J + O ((H tj y) . (12) 
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The leading term in eq. QT2D arises by using the diagonal part of the Hamiltonian 
in eq. ([U])- The off-diagonal part of the Hamiltonian yields sub-leading terms. The 
resulting equation suggests a damping of the non-diagonal elements which decrease 
with increasing I. 



Let us first take the original Hamiltonian with ^-dependent couplings and study 
its evolution with the Hamiltonian flow equations. We improve the scheme in the 
second half of the paper by adding a term which starts out as zero for I = but 
grows during the evolution. The generalized Hamiltonian T~i(£), together with the 
initial conditions, has the following form: 

H(£) = t{l)J z + V{l){Jl + Jl), 
£(0) = Co, 

V(0) = V . (13) 
Defining rj as r\ = [TCdjH] with TCd = £,{£)J Z we obtain: 

77 = 2£V(Jl -J 2 _). (14) 
By inserting this r\ into the flow equation we find: 
dH 



= 4ttV 2 (2J I (T 1 -Jl)-J;)-4eV(Jl + J' i _). (15) 

The right hand side of this equation contains an hermitian operator containing three 
quasi spins oc , i.e. a term different from the original evolving Hamiltonian (13). 
This new term is approximated by linearizing around the Hartree-Fock operator 
expectation values: 

AB = (A — (A))(B - (B)) - (A)(B) + (A)B + A(B) (16) 

which yields: 

J 3 Z ^3(J Z ) 2 J Z . (17) 

A linearization of J| is certainly meaningful when the fluctuations of J z around its 
expectation value are small, i. e. for large quasi spins. This is the case for the lowest 
lying levels, if the particle number is large. 

Differentiation of the ansatz for TC(£) in eq. (|1^) gives the left side of the evolution 
equation 

-J z + Zr(J + + J-) (18) 



dl d£ 

which we compare with the approximate result for the commutator from eq. 

We see that the ansatz for Tt{£) in eq. fllBp together with the approximation for 
the cubic operator product indeed leads to a closed system of first order differential 
equations for the coefficients in front of the quasi-spin operators. 

H = -4£V 2 LV(iV-l) + l], 

% = ~^ 2 V- (19) 
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The second equation implies that the magnitude of the off-diagonal matrix element 
V decreases in the course of the evolution, independently of its sign. Combining the 
two equations one readily obtains an invariant of the evolution: 

e(£) - [N(N - 1) + 1]V 2 (£) = e - [N(N - 1) + IK 2 . (20) 

Using V(oo) = one finds the limiting value of £ and the splitting A between the 
ground state and first excited state : 



f(oo) = ^oVT^ 



N(N-1) + A 

A = (l|W(oo)|l)-(0|W(oo)|0) 

= £(oo). (21) 

In the large iV-limit with finite interaction strength the splitting A between the 
ground state and first excited state has the form of the RPA approximation H 



2 



A rp a = £oVl - e 

e = iV/. (22) 

In nuclear physics one studies systems with intermediate particle number. Nei- 
ther few nor many body methods are directly applicable. In nuclei the relevant 
number of particles is given by the number of particles outside of closed shells, 
which is iV < 2 J + 1, where J is the angular momentum of the valence shell. We 
want to obtain a reliable result also for small values of N. To this end we propose to 
add a new operator to the original Hamiltonian, which improves the accuracy of the 
evolution equations. The new term oc has been chosen to have the same symme- 
try as the original Hamiltonian under reflections. From the previous flow equations 
such an improvement is suggested, since a term oc J| was generated on the right 
hand side of the previous evolution equation (|To]). Adding new operators generated 
from the original Hamiltonian by the commutators is an established way to improve 
the method. Obviously symmetry provides an important guideline for the choice of 
new operators to be included in the Hamiltonian 7i{€). Usually the commutators 
generate additional terms which must be eliminated, since one can obtain a closed 
set of differential equations only if every operator on the right side of the evolution 
equation has a counterpart on the left-hand side of eq. @. We choose: 

W) = + V{l){J\ + Jl)+i(£)J 3 z (23) 

with the initial conditions 

m = £o; V(0) = V - 7 (0) = 0. (24) 

Proceeding along similar lines as before we define rj as 

V=IZJ Z ,H] (25) 
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slightly deviating from the commutator of [Hdj'H]- Otherwise terms with higher 
powers of quasi spins would be produced. The explicit calculation of rj leads to the 
same anti-hermitian operator product as before. 

r) = 2{;v(Jl- J -)- ( 26 ) 

Inserting this rj into the flow equation we get 



^ = A^V 2 (2J Z (J 2 - J 2 ) - J z ) - 4£V(4 + Jl) + 2tV 1 [(J 2 + - J 2 ), Jl\. (27) 

We see a new operator product appearing from the commutator of i] with the 7 
term in Ti. Using the reduction scheme for products of operators described above, 
we find 

[(Jl-J 2 _),J!]^-6(J z nJl + J 2 _) (28) 

This way we regain a closed system of first order evolution equations for the coeffi- 
cients in front of the quasi-spin operators 



di 
dV 

~dJ 
~di 



A^V 2 [2J(J + 1)-1], 
-4£V - 12^ 7 (J,) 2 , 

-8£V 2 . (29) 



We introduce the abbreviations 

r = 2J 2 - 1 = 2J(J+ 1) - 1, 

s = \{J Z )\ 

t = ^. (30) 
r 

Combining the first and third evolution equation we find the invariant 

£ + 2 = ^, (31) 
r 2 r 

where we employ the initial conditions on the right hand side. By using this relation 
we eliminate 7 as independent variable. A second more complicated invariant arises 
from all three equations and has the form 

? + rV l _ t ^ Q _ £)2 = £2 + rV 2. (32) 

The allowed values of ^ and V lie on a hyperbola as shown in fig. 1 . For the plot we 
choose e = 0.8 and t — 3, which is approximately true for large N. At the beginning 
of the evolution at £ = we start on the hyperbola at £ = £ an d V = V . Since the 
gradient d^jdt > the magnitude of £ is growing and V decreases. Thus the point 
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Figure 1: Flow in the £ — V plane as denned by eqs. (29) and (32) 



V = 0, where the right hand sides of all three evolution equations vanish, is a stable 
fix point. The variables £ and 7 converge towards the fix point values 



£(00) 



7(00 J 



t- 1 

2£o 
r(t-l) 



'1-3 



t- 1 

2t 



'1-3 



t - 1 
2t 



(33) 
(34) 



The above expressions depend on the particle number N through the constants r, s, t 
in a rather complicated way. Due to the presence of the J| term in the Hamiltonian 
the splitting A between the ground state and first excited state has the modified 
form: 



(l|W(oo)|l)-<0|W(oo)|0) 



f(oo) + 7(oo)( 



3N 2 



(35) 



In fig. 2 we show the energy gap A for three different interaction strengths e = 
0.4, 0.6, 0.8. The curves with stars are the result of an exact diagonalization of the 
Hamiltonian. The full drawn curves give the improved result from the evolution 
equations (^). For comparison we also show the result of the evolution equations 
(P^|) with long dashes and the RPA result ( |22"D with short dashes. One sees that 
the method is very accurate for weak interactions, where perturbation theory also 
does a rather good job. In the case of strong coupling there are still some deviations 
for intermediate particle numbers. The larger the coupling the closer one is to 
the transition to a deformed ground state, therefore also the fluctuations become 
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stronger and our expansion of the new operators in powers of fluctuating operators 
is not so good anymore. All three calculations approach the RPA result for large 
particle numbers shown with short dashes. In the limit of large N one can also 
demonstrate analytically that the improved method converges to the RPA result. 

Applying the evolution equations to the Lipkin model one encounters the problem 
that ever new operators are generated in the course of setting up the evolution 
equations. This is a general feature of spin models. Similar features occur in the 
spin boson problem. The original Lipkin Hamiltonian has a rather simple tridiagonal 
structure, which is lost during the evolution for finite £, if one does not enforce it by 
suitable approximations. Recently Mielke [ 10 1 has constructed an 77-operator which 



preserves the simple structure of the Lipkin Hamiltonian without approximations, 
but still reduces the magnitude of the off-diagonal matrix elements. This method is 
also applicable beyond the stability of the Hartree-Fock state (e > 1). 

In conclusion, we have set up Hamiltonian flow equations for the Lipkin model. 
By infinitesimal unitary transformation we have diagonalized the Hamiltonian ap- 
proximately making reasonable truncations for the new operators. We have pre- 
sented an improved scheme which successfully interpolates between perturbation 
theory for small iV and RPA for large N. The method of evolution equations has 
been shown to be a promising new tool also for nuclear physics. It opens up new 
perspectives for approximations in many body problems which still have to be ex- 
plored in greater detail. Thus a systematic approach to select the relevant set of 
operators for the evolution equations can be established. 
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Figure 2: Energy gap A between ground state and first excited state for various 
strengths e = 0.4, 0.6, 0.8. Stars denote the exact result, solid line shows the im- 
proved result from the evolution equations, while the long-dashed line gives the 
result of the linearized scheme. The short-dashed line is the RPA result. 
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